Parametrically controlling solitary wave dynamics in modified Kortweg-de Vries 
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We demonstrate the control of solitary wave dynamics of modified Kortweg-de Vries (MKdV) 
equation through the temporal variations of the distributed coefficients. This is explicated through 
exact cnoidal wave and localized soliton solutions of the MKdV equation with variable coefficients. 
The solitons can be accelerated and their propagation can be manipulated by suitable variations of 
the above parameters. In sharp contrast with nonlinear Schrodinger equation, the soliton amplitude 
and widths are time independent. 
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Modified Kortweg-de Vries (MKdV) equation mani- 
fests in diverse areas of physics 0, 0, 0, 0, IE @ F° r 
example, it appears in the context of, electromagnetic 
waves in size-quantized films, van Alfven waves in colli- 
sionless plasma , phonons in anharmonic lattice 8] , in- 
terfacial waves in two layer liquid with gradually varying 
depth 0, transmission lines in Schottky barrier [lCf. ion 
acoustic solitons ITllJl3.IT^ | , elastic media 0] , and traf- 
fic flow problems [la, [l6| • It is an integrable dynamical 
system with an infinite number of conserved quantities; 
the solutions of this equation are well studied [13, EH ■ 

Recently non-linear equations with variable coefficients 
have attracted considerable attention in the literature. 
Nonlinear Schrodinger equation (NLSE) with variable 
non-linearity and dispersion is relevant to bot h op tical 
fibers and Bose-Einstein condensates 0, |2(], l2ll |22| . 
Nonlinear Schrodinger equation with source, having dis- 
tributed coefficients like variable dispersion, variable 
Kerr nonlinearity and gain or loss .is applicable to asym- 
metric twin-core optical fibers 23,24]. It has been shown 
that, solitons can be compressed and their dynamics ef- 
fectively controlled through these variable parameters. 
The Kort weg -de Vries (KdV) equation with variable co- 
efficients |25| has been studied recently in the context of 
ocean waves, where the spatio-temporal variability of the 
coefficients are due to the changes in the water depth and 
other physical conditions. The fact that, MKdV equa- 
tion is relevant to hydrodynamics and a variety of phys- 
ical phenomena, it is natural to expect the possibility of 
temporal variations in the equation parameters occurring 
in the same. Furthermore, for propagating solitons, the 
first integral of the MKdV equation yields NLSE with a 
source, making it imperative to investigate the effect of 
the temporal variation of the distributed parameters on 
the solitary wave solutions of this dynamical system. 

The goal of the present paper is to study the effect of 



the variable coefficients on the solution space of MKdV 
equation, both for positive and negative cases. We find 
that the effect of distributed coefficients on the soliton 
dynamics of MKdV equation is quite different than that 
of NLSE. In case of NLSE, the amplitude and width are 
affected by the time dependence of the distributed coef- 
ficients. This leads to compression of solitary waves in 
NLSE. Through explicit construction, it is shown that, 
solitary waves in this system can be effectively controlled 
through the equation parameters. The solitons can be ac- 
celerated and manipulated by suitable variations of the 
above parameters. However the width and amplitudes 
are not amenable for manipulation and control. 

We consider the modified KdV equation with variable 
coefficients in the form 



u t + a(t)u x - I3(t)u 2 u x 4- j(t)u xxx = 0, 



(1) 



where j(t), a(t) and (3(t) are time dependent variables. 
Although, the first derivative term in the field variable 
can be removed by a suitable change of the coordinate 
frame, the same has been kept here explicitly to contrast 
its effect with the nonlinear and dispersion terms. We 
consider the ansatz solution of the form, 



u = A 1 (t)g[oj(x,t)]+A (t), 



(2) 



where, ui(x,t) = f(t)x — h(t). The variable coefficient 
MKdV equation can be mapped to Jacobi elliptic equa- 
tion: 



g = Pg + 2Qg\ 
having a conserved quantity, 

(g) 2 = Pg 2 + Qg 4 



(3) 



(4) 



Here prime indicates differentiation with respect to the 
argument to, and P and Q are constants. Substituting 
the ansatz in Eq.Q and using the relations and l@J, 
we collect the coefficients of g a and g a g (where a=0,l,2) 
to find the consistency conditions: 
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5° : d t A = 0, 



(5) 
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9 ■ d t Ax = 0, 



(6) 



Axdtuj+A^dxU+A^it^dM^-AlA^dxUj = 0, 
gg : -2A Q A\p{t)d x u = 0, 



Below we analyze some explicit solutions and corre- 
sponding parameter ranges. For the sake of specificity 
we consider (3{t) > 0. 

Case I: With g = cn(uj(x, t)) one finds, the cnoidal wave 
solution as, 



(7) 
(8) 



g 2 g : -A\(3{t)d x w + ^Q^^tof = 0, (9) 

For obtaining solutions to the above set of equations, we 
require, Aq = 0, and Ai =constant. From Eq.@, further 
simplification yields, 



fit) 



Aim 



(10) 



6Q7(*)' 
and from Eq.J7Jl we get, 

O t h(t) = xd t f(t) + a(t)f(t) + Pj(t)f(t). (11) 

Condition requires that, f(t) should be a constant 
so that the term containing x vanishes. This implies 
f3(t)/j(t) — k, where k is a constant. We then have the 
relations, 




and, 



6Q 



(12) 



(13) 



Hence, the exact travelling wave solution can be writ- 
ten in the form, 



u = A x g 



1 A\k 



It is worth noting that, A\ is unconstrained and controls 
the width of the solution. Unlike the case of NLSE the 
amplitude and width are independent of time. Since, the 
solution involves k, positive and negative MKdV equa- 
tions have different type of solutions, g can be any of the 
twelve Jacobi elliptic functi ons, with the modulus param- 
eter to 2 (0 < to 2 < 1) [2l,l23- The following arc some 
identities of the Jacobi elliptic functions, which are used: 
cn 2 (w, m) + sn 2 (w, m) = 1, 
dn 2 (w, to) + m 2 sn 2 {w, to) = 1, 
sn (w, to) = cn(w, m)-dn(w, to), 
an (ui,to) = —sn(w,m)-dn(w,m), 
dn (w,m) — —m 2 sn(w,m)-cn(w,m). 
For to = 1, 

cn(w, 1) = dn(w, 1) = sech(w) and sn(w, 1) = tanh(w). 



u = A\cn 



A\k 
—6m 2 



x- [a{t) + 



7(*)(2r 



1)A 2 K. 



-6m 2 



dt 



(15) 

where, j(t) < 0, P = (2m 2 - 1) > 0, Q = -m 2 < 
and to 2 > 1/2. This corresponds to the positive MKdV 
equation. In the case, when m 2 = 1/2, P = 0, and j(t) 
does not affect the solution. For, to 2 = 1 case, we have 
the exact solitary wave solution of the form, 



Ai sech 



>-A 2 n 



Ht) 



l(t)A 2 lK 
6 



]dt 



(16) 

Fig. J3> depicts the temporal evolution of the above 
bell shaped localized solution. For illustrative purpose we 
have considered two different cases where, a(t) = and 
a(t) = —3t 3 cos(t 3 ). Fig. (0) depicts the same solution 
when j(t) and a(t) have polynomial time dependence. 
One clearly sees that the temporal variations of 7 and a 
can effectively modulate and control the propagation of 
the solitons. 
Case II: g = sn(w, to) 

We now study the cases where g = sn(uj(x,t)), for 
which the solution corresponds to negative MKdV equa- 
tion: 



u = A\sn 



IA(k 
6to 2 



7 (t)(r 



; ..-u,<> t1)^k u 
[a{t) — — \dt 



6m 2 



(17) 

Here j(t) > 0, Q = to 2 > and P = -(to 2 + 1) < 0. 
For to 2 = 1, we have the kink type solitary wave solution 



u = Aitanh 



*A\k 
~6~ 



a(t) - lMH ]dt 



(18) 

Fig. © depicts the kink solution in the presence of 
time dependent dispersion. The soliton motion can be 
controlled through the external parameters. 
Case III: g = dn(w, m) 

We get exact solitary wave solution of positive MKdV 
equation, in the case when g = dn(w(x,t)): 



u = A\dn 



—A\k 



. 7 (^)(2-to 2 )A 2 K i ^ 
a(t) \dt 

(19) 
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Here, j(t) < and P = (2 - to 2 ) > 0. For, m 2 = 1, 
we get bell shaped solitary wave solution as, 



u = A\dn 



-A\k 



G 



a{t) \dt 



(20) 

In conclusion, MKdV equation with time varying co- 
efficients has solitary waves solutions, provided the tem- 



poral variations of the coefficients are of the form given 
in the text. The temporal variation of these parameters 
allow effective control of the solitary wave profile. These 
continuous waves and localized solutions can be made to 
accelerate. The amplitude and widths are not modulated 
by the distributed coefficients. The induction of time de- 
pendent u x term allows us to control the motion of the 
solitons more efficiently. 
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FIG. 1: Propagating localized solitary wave solution of MKdV equation with g = cn(x,i), where 7 = cos(f), k = —24, A\ = 1, 
P = 1, Q = -1 and m = 1; left a(i) = and right a(t) = -3t 3 cos(t 3 ). 




FIG. 2: Propagating localized solitary wave solutions for g = cn(x,t), where j(t) = 3£ 2 , k = — 24,^4 1 = 1, P = 1, Q = —1 and 
m = 1; left a(i) = and right a(t) = -t 9 . 
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FIG. 3: Kink type solitary wave solution of the MKdV equa- 
tion, where j(t) = cos(£), a(t) = 5t 4 , k — —48, A\ = 1, 
P = 1, Q = — 1 and m = 1. 



